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Abstract 

The  specific  character  of  biological  enzyme  catalysts  enables  combined  fuel  and  oxidant  channels  and  simplified  non-compartmentalized  fuel 
cell  assemblies.  In  this  work,  a  micro  structured  enzymatic  biofuel  cell  architecture  is  proposed,  and  species  transport  phenomena  combined  with 
consecutive  chemical  reactions  are  studied  computationally  in  order  to  provide  guidelines  for  optimization.  This  is  the  first  computational  study  of 
this  technology,  and  a  2D  CFD  model  for  species  transport  coupled  with  laminar  fluid  flow  and  Michaelis-Menten  enzyme  kinetics  is  established.  It 
is  shown  that  the  system  is  reaction  rate  limited,  indicating  that  enzyme  specific  turnover  numbers  are  key  parameters  for  biofuel  cell  performance. 
Separated  and  mixed  enzyme  patterns  in  different  proportions  are  analyzed  for  various  Peclet  numbers.  High  fuel  utilization  is  achieved  in  the 
diffusion  dominated  and  mixed  species  transport  regimes  with  separated  enzymes  arranged  in  relation  to  individual  turnover  rates.  However,  the 
Peclet  number  has  to  be  above  a  certain  threshold  value  to  obtain  satisfying  current  densities.  The  mixed  transport  regime  is  particularly  attractive 
while  current  densities  are  maintained  close  to  maximum  levels.  Optimum  performance  is  achieved  by  mixed  enzyme  patterning  tailored  with 
respect  to  individual  turnover  rates,  enabling  high  current  densities  combined  with  nearly  complete  fuel  utilization. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

In  recent  years,  the  demand  on  the  power  supplies  of  portable 
electronic  equipment  such  as  laptops  and  cellular  phones  has 
substantially  increased  due  to  additional  technical  features  and 
miniaturization  efforts.  Well-developed  lithium  ion  battery  tech¬ 
nology  currently  powers  most  of  these  applications.  It  is  unlikely, 
however,  that  improvements  in  battery  technology  will  keep  pace 
with  the  accelerating  power  demand  [1].  In  order  to  enable  fur¬ 
ther  advances,  other  alternatives  such  as  fuel  cells  have  to  be 
considered.  The  high  cost  associated  with  the  development  of 
new  small-scale  fuel  cell  power  supplies  are  surpassed  by  the 
growing  portable  power  needs  of  the  mobile  workforce  and  the 
major  industries  that  now  support  them  [2] .  Direct  methanol  fuel 
cells  exhibit  relatively  high  power  and  energy  densities  [3],  but 
there  are  several  problems  associated  with  these,  for  instance 
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methanol  crossover  and  slow  anode  kinetics  reduce  efficiency 
and  open-circuit  voltage  [4-5].  Recently,  this  technology  has 
been  challenged  by  emerging  microscale  biofuel  cell  designs. 

A  biofuel  cell  is  not  necessarily  a  fuel  cell  that  uses  a  biolog¬ 
ical  fuel,  although  that  might  very  well  be  the  case.  Biofuel  cells 
always  involve  biological  molecules  or  living  cells  to  catalyze 
the  chemical  reactions,  replacing  conventional  transition  metal 
catalysts  such  as  platinum.  The  underlying  concept  is  to  control 
electrochemical  processes  occurring  in  nature  to  provide  useful 
electrical  current.  Chemical  energy  can  be  converted  to  electri¬ 
cal  energy  by  coupling  the  biocatalytic  oxidation  of  organic  or 
inorganic  compounds  to  the  chemical  reduction  of  an  oxidant 
in  a  controlled  environment.  There  are  two  main  types  of  bio¬ 
fuel  cells;  microbial  biofuel  cells  and  enzymatic  biofuel  cells. 
Microbial  fuel  cells  utilize  entire  living  cells  or  microorganisms 
combined  with  redox  intermediates  to  catalyze  the  oxidation  of 
a  fuel.  Microbe  catalyzed  processes  involve  several  fairly  com¬ 
plex  steps  with  various  intermediates  [6]  and  are  quite  inefficient 
[7].  In  an  enzymatic  biofuel  cell,  the  chemical  reactions  are  cat¬ 
alyzed  by  biological  redox  enzymes  that  can  be  separated  and 
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Nomenclature 

Ac 

cross-sectional  area  (m2) 

c 

concentration  (molm-3) 

D 

solute  diffusivity  (m2  s-1) 

Dh 

hydraulic  diameter  (m) 

Da 

Damkohler  number 

F 

Faraday’s  constant  (C  mol-1) 

FU 

fuel  utilization 

8 

constant  of  gravity  (m  s-2) 

h 

channel  height  (m) 

i 

current  density  (A  m-2) 

J 

species  molar  flux  (mol  s-1  m-2) 

k 

rate  constant  (s-1) 

^cat 

turnover  number  (s-1) 

Km 

sensitivity  (molm-3) 

L 

axial  length  of  electrode  (m) 

Lc 

characteristic  length  (m) 

M 

molar  mass  (kg  mol-1) 

Na 

Avogadro’s  number  (mol-1) 

P 

pressure  (Pa) 

Pe 

Peclet  number 

R 

turnover  rate  (mol  s-1  m  2) 

Re 

Reynolds  number 

t 

time  (s) 

T 

temperature  (K) 

u 

bulk  velocity  (ms-1) 

u 

mean  bulk  velocity  (ms-1) 

Vout 

output  voltage  (V) 

V 

volumetric  flow  rate  (m3  s-1) 

V 

specific  activity  (mol  kg-1  s-1) 

w 

channel  width  (m) 

Wout 

output  power  (W) 

ftp 

pumping  power  (W) 

X 

direction  of  main  flow  (axial) 

y 

direction  across  channel  (transverse) 

z 

number  of  electrons 

Greek  letters 

fi 

kinematic  viscosity  (kg m  ]s  l) 

P 

density  (kgm-3) 

Subscripts 

in 

inlet  value 

max 

maximum  value 

out 

outlet  value 

ref 

reference  value 

0 

bulk  value 

Superscripts 

A 

methanol  dehydrogenase  (MDH) 

ADH 

alcohol  dehydrogenase 

B 

formaldehyde  dehydrogenase  (FalDH) 

C 

formate  dehydrogenase  (FDH) 

purified  from  suitable  organisms,  thus  extracting  the  actual  part 
of  the  cell  that  enables  catalytic  operation  [8].  Enzyme  catalysts 
are  specific  to  particular  substances  (fuels)  and  the  presence  of 
other  substances  does  not,  in  general,  impact  the  rate  of  cataly¬ 
sis.  This  specificity  of  enzyme  catalysis  enables  the  combination 
of  fuel  and  oxidant  streams  in  a  single  manifold  [8],  with  multi¬ 
ple  benefits  in  the  context  of  fuel  cell  design  and  operation:  no 
proton  exchange  membrane  is  needed,  sealing,  manifolding  and 
fluid  delivery  requirements  are  greatly  reduced,  ohmic  losses 
and  water  management  issues  associated  with  proton  exchange 
membrane  fuel  cell  (PEMFC)  technologies  are  eliminated  [9], 
as  are  two-phase  transport  issues  associated  with  gas  diffusion 
layers  [10].  With  fuel  and  oxidant  streams  fully  mixed  in  the 
same  channel,  fuel  crossover  is  also  not  a  problem.  Cost  reduc¬ 
tion  is  possibly  the  main  aspect  that  can  make  biological  fuel 
cells  competitive  for  small-scale  power  supplies.  Firstly,  the  pro¬ 
ton  exchange  membrane  is  a  costly  component  of  conventional 
fuel  cell  assemblies  [11].  Secondly,  the  cost  of  platinum,  which 
is  already  high,  is  expected  to  increase  further  with  increased 
demand  and  production  [12].  Biological  enzymes,  on  the  other 
hand,  are  derived  from  naturally  abundant  organisms  and  can 
be  produced  using  low-cost  fermentation  techniques  [13].  Other 
benefits  are  operation  near  or  at  room  temperature  and  compact 
units  produced  by  inexpensive  and  well-established  microfab¬ 
rication  methods.  Channels  down  to  10  p,m  in  diameter  can  be 
fabricated,  greatly  improving  the  active  surface-to- volume  ratio 
of  the  device  and  consequently  its  power  density,  which  remains 
quite  limited  with  current  biofuel  cell  designs.  As  with  all  fuel 
cell  technologies,  however,  stability  is  a  key  issue  for  biofuel 
cells  [13].  Enzymatic  biofuel  cells  must  withstand  extended  use, 
varied  duty  cycles,  a  range  of  environmental  conditions  and  must 
reach  lifetimes  satisfying  the  demands  of  advanced  applications, 
preferably  on  the  order  of  years.  Possible  applications  range 
from  miniaturized  portable  electronic  equipment  [2]  and  sen¬ 
sors  to  integrated  lab-on-chip  power  supplies  [14]  and  advanced 
in  vivo  diagnostic  medical  devices  with  reactants  available  in 
the  ambient  environment  [13]. 

In  most  enzymatic  biofuel  cells  developed  to  date,  enzymes 
are  placed  in  solution  and  the  electron  transfer  from  the  enzyme 
biocatalysts  to  the  electrode  is  mediated  by  a  redox  couple.  The 
rate  of  electron  transfer  is  generally  limited  by  the  ability  of 
the  redox  species  to  diffuse  back  and  forth.  Several  methods  to 
improve  electron  transfer  rates  by  immobilization  of  enzymes 
on  electrode  surfaces  have  been  presented,  for  example:  cova¬ 
lent  enzyme-electrode  polymer  linkages,  layered  electrodes 
obtained  by  functional  monolayers  and  non-covalent  coupling 
by  hydrophobic/hydrophillic  or  affinity  interactions  [15].  Elec¬ 
tron  mediators  (co-factors)  extracting  electrons  from  the  active 
sites  of  the  enzymes  are  necessary  for  efficient  flow  of  electrical 
current  [8].  Alcohol  dehydrogenase  enzymes  were  immobilized 
in  Nation  by  Moore  et  al.  [14]  and  their  ethanol  fuelled  anode 
reached  an  open-circuit  voltage  of  0.34  V  and  current  densities 
of  53  jxA  cm-2  using  an  external  platinum  cathode.  The  highest 
open-circuit  voltage  (0.7 1  V)  and  power  density  (1.5  mW  cm-2) 
for  a  methanol/oxygen  biofuel  cell  was  observed  when  the  bio¬ 
catalyst  was  immobilized  close  to  the  anode  and  utilized  its 
electron  mediator  [16].  To  achieve  this,  complete  oxidation  of 
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methanol  to  carbon  dioxide  by  three  different  enzymes  was 
engaged.  Patterning  multiple  enzyme  electrodes  to  harness  con¬ 
secutive  reactions  represents  a  significant  opportunity  for  biofuel 
cell  technologies,  particularly  with  respect  to  fuel  utilization. 

The  biocatalysts  of  interest  here  are  covalently  linked  directly 
to  the  electrode  surface  via  electron  mediators  (cofactor)  that  are 
immobilized  in  contact  with  the  active  enzyme  sites.  Biomimetic 
dyes  are  used  as  electron  mediators;  they  insert  into  the  binding 
site  on  the  enzyme  and  bind  specifically  and  irreversibly  without 
interfering  with  the  catalytic  activity  of  the  enzyme.  An  illustra¬ 
tion  of  the  proof-of-concept  microfluidic  biofuel  cell  is  shown 
in  Fig.  1,  with  a  schematic  depiction  of  the  surface-bound  bio¬ 
catalysts.  The  simple  flat-plate  microchannel  design  facilitates  a 
high  surface-to-volume  ratio,  effectively  one-dimensional  flow 
behavior,  rapid  cross-stream  diffusion  and  straightforward  sur¬ 
face  preparation.  Prior  to  assembly,  an  electrically  conducting 
material  is  deposited  on  the  interior  surfaces  to  form  the  anode 
and  cathode,  and  the  enzymes  are  tethered  directly  to  these  lay¬ 
ers. 

The  focus  of  this  study  is  to  model  species  transfer  associ¬ 
ated  with  heterogeneous  chemical  reactions  and  enzyme  kinetics 
based  on  the  microchannel  geometry  in  Fig.  1.  The  first  aim  is 
to  determine  whether  the  process  is  diffusion  limited  or  reac¬ 
tion  rate  limited.  Secondly,  various  enzyme-electrode  patterns 
coupled  with  coherent  bulk  velocities  are  investigated  in  order 
to  realize  efficient  fuel  cell  operation.  This  is  the  first  computa¬ 
tional  study  of  this  technology,  and  the  results  are  intended  to 
provide  guidelines  for  the  design  and  fabrication  of  microfluidic 
biofuel  cells  exploiting  consecutive  reactions.  The  long-term 
objective  is  to  integrate  the  forefront  of  molecular  biology  with 
advanced  microfabrication  techniques  to  produce  practical  bio¬ 
fuel  cell  systems  with  novel  architectures  that  result  in  advanced 
functionality. 

2.  Method 

2.7.  Enzyme  kinetics 

An  enzyme  converts  one  chemical  species  (the  substrate) 
into  another  one  (the  product).  The  substrate  concentration  dis¬ 
tribution  in  the  combined  fuel  and  oxidant  feed  stream  of  a 


microfluidic  biofuel  cell  depends  on  the  rate  of  the  heteroge¬ 
neous  chemical  reactions  taking  place  on  the  biocatalytic  elec¬ 
trode  surfaces.  For  now,  we  ignore  the  coupled  effects  of  oxidant 
transport  and  cathode  kinetics  and  instead  focus  on  the  anodic 
process.  The  Damkohler  number  represents  the  ratio  of  reaction 
velocity  to  diffusive  velocity  defined  as 


In  this  equation,  k  is  a  rate  constant  for  heterogeneous  reactions, 
co  the  bulk  concentration,  v  the  overall  order  of  reaction,  D  the 
diffusion  coefficient  and  Lc  is  a  characteristic  length.  There  are 
two  common  scenarios: 

(1)  Da  1:  diffusion  limited  scenario 

The  reaction  rate  is  fast  compared  to  substrate  diffusion 
to  the  catalyst  site.  Substrate  transport  is  occurring  at  its 
maximum  rate  with  approximately  zero  concentration  at  the 
reaction  surface. 

(2)  Da<^  1 :  reaction  rate  limited  scenario 

Diffusion  of  species  is  fast  compared  to  the  reaction  rate. 
The  concentration  of  substrate  at  the  reaction  surface  is 
approximately  equal  to  the  bulk  concentration. 

To  determine  if  the  biofuel  cell  is  diffusion  limited  or  reaction 
rate  limited,  an  appropriate  system  to  model  is  the  well-known 
horse  liver  alcohol  dehydrogenase  (ADH)  enzyme  kinetics  for 
ethanol  oxidation  to  acetaldehyde.  The  ADH  catalyzed  reaction 
is  a  one-step  process,  namely 

CH3CH2OH^*CH3CHO  +  2H+  +  2e“ 

Wild  type  horse  liver  ADH  has  a  maximum  turnover  number  for 
catalytic  ethanol  oxidation  of  k^u  =  155  s-1  for  pH  from  8 
to  10  [17].  It  is  a  dimer  with  an  irregular  shape  which  makes 
it  difficult  to  measure  its  size.  The  unit  cell  of  a  P2i  space 
group  of  a  crystal  structured  form  of  ADH  is  reported  to  be 
55.05  x  73.63  x  92.49  A  (0=  102.6°)  [18].  As  a  rule-of-thumb, 
its  molecular  cross-sectional  area  can  be  estimated  to  be  around 
Ac  =  10“ 16  m2.  If  we  assume  that  the  electrode  surface  cover¬ 
age  of  enzymes  is  complete,  the  maximum  turnover  rate  for  the 
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heterogeneous  reaction  can  be  calculated  as 

flmax  =  *cat[E]  =  (2) 

AcA^a 

The  surface  concentration  of  enzymes  is  in  our  case  equal  to 
[E]  =  (A,JVa)-1  •  Thus,  we  get 

dADH  ^ca?H  _ _ 

max  AcNa  10-16  x  6.02  x  1023 

=  2.6  x  10-6molm-2  s_1 


The  one-dimensional  diffusive  species  flux  is  defined  as 


dc 

J  =  -D— 
dx 


(3) 


If  we  consider  a  microchannel  h-  20  |mm  in  height  with  a  bulk 
ethanol  concentration  of  500  mM,  the  maximum  diffusive  flux 
of  ethanol  becomes 


J  =  -D—  =  10“9 
h 


0.5  x  103 
20  x  10"6 


=  2.5  x  10-2  mol  m-2  s_1 


Thus,  species  diffusion  is  about  four  orders  of  magnitude  faster 
than  the  reaction  rate,  indicating  that  the  system  is  reaction  rate 
limited  ( Da< ^  1).  Even  if  the  turnover  rate  was  improved  by 
a  factor  of  100,  the  anodic  process  would  still  be  reaction  rate 
limited. 

The  methanol  fuelled  biofuel  cell  oxidizes  methanol  to  car¬ 
bon  dioxide  via  the  simplified  three-step  process 


CH3OH^>HCHO+2H+  +  2e~ 


information  available  in  enzyme  databases.  The  relatively  low 
MDH  activity  is  strongly  enhanced  by  a  50  kDa  activator  protein 
from  the  same  organism,  resulting  in  a  40-fold  increase  in  the 
MDH  turnover  rate  (included  in  the  tabulated  value)  [20] .  A  pro¬ 
tein  purified  from  Methylococcus  capsulatus  has  been  shown  to 
improve  the  specific  activity  of  FalDH  [23] .  There  are  also  higher 
specific  activities  available  for  FDH  purified  from  Pseudomonas 
oxalaticus  with  cofactor  methylene  blue  or  NAD  [24].  Turnover 
numbers  are  normally  not  available  in  the  literature,  but  can  be 
calculated  based  on  the  specific  activity,  v,  as  described  by  the 
Michaelis-Menten  model  for  first-order,  steady-state  enzyme 
kinetics  [25], 


PmaxQ) 
A  m  +  Cq 


(4) 


The  sensitivity  Km  is  a  measure  of  the  changes  in  reaction 
velocity  when  the  substrate  availability  is  reduced,  in  fact,  Km 
represents  the  substrate  concentration  value  at  which  the  reac¬ 
tion  velocity  is  half-maximum,  v  =  vmax/2.  Catalytic  oxidation 
of  primary  alcohols  by  alcohol  dehydrogenases  using  cofac¬ 
tors  in  fact  resembles  two- substrate  reactions  described  by  the 
Theorell-Chance  mechanism,  where  a  ternary  complex  is  pre¬ 
sumably  formed  but  its  breakdown  to  yield  the  first  product  is 
very  fast  so  that  the  ternary  complex  is  kinetically  insignif¬ 
icant  [26].  With  the  cofactor  bound  to  the  enzyme,  one  can 
assume  that  the  binding  of  substrate  is  the  rate-limiting  step  as 
in  an  ordinary  one- substrate  mechanism,  expressed  in  the  clas¬ 
sic  Michaelis-Menten  model  of  two-stage  enzyme-substrate 
action, 


HCHO  +  H20F^5IHC00“  +  3H+  +  2e“ 


e  +  s24esAe  +  p 

k-\ 


HCOO_™C02  +  H+  +  2e“ 

The  turnover  numbers  and  average  sizes  of  the  three  cat¬ 
alytic  enzymes  participating  in  this  process  are  expected  to  be 
of  the  same  order  of  magnitude  as  for  ADH,  and  all  systems 
are  therefore  expected  to  be  reaction  rate  limited.  These  new 
enzymes  appear  less  frequently  in  the  literature  and  there  is  rel¬ 
atively  little  appropriate  data.  Table  1  presents  a  summary  of  the 


The  single  substrate  S  binds  reversibly  to  the  enzyme  E  in  the 
first  stage,  and  the  single  product  P  is  formed  in  the  second 
stage.  The  sensitivity  in  Eq.  (4)  combines  the  effects  of  forward 
rate  constant  k\  ,  back  rate  constant  k-\  and  product  formation 
rate  constant  k 2.  Fig.  2  shows  a  typical  Michaelis-Menten  curve 
with  normalized  variables  v/vmax  and  co/Km.  If  the  sensitivity 
is  small  relative  to  the  bulk  concentration,  co,  the  specific  activ¬ 
ity  is  close  to  its  maximum  value.  This  implies  that  fuel  cell 


Table  1 

Enzyme  data  valid  for  enzymes  in  solution 


MDH  [20] 

FalDH  [21] 

FDH  [22] 

ADH  [17-19] 

EC  number 

1.1.1.244 

1.2.1.46 

1.2.1. 2 

1.1. 1.1 

Organism 

B.  methanolicus 

P.  putida 

P.  putida  sp.  101 

E.  caballus 

Substrate 

Methanol 

Formaldehyde 

Formate 

Ethanol 

Product 

Formaldehyde 

Formate 

C02 

Acetaldehyde 

Cofactor 

NAD 

NAD 

NAD 

NAD 

Number  of  subunits 

10 

4 

6 

2 

Molecular  weight  (kDa) 

430 

250 

280 

- 

Specific  activity  (ixmolmin-1  mg-1) 

10.4 

8.3 

2.43 

- 

Max  substrate  cone.  (M) 

0.5 

- 

- 

- 

Required  ions 

Mg2+,  Zn2+ 

- 

- 

- 

Optimum  temperature  (°C) 

- 

30-40 

40 

- 

Optimum  pH 

- 

8.5 

8 

8-10 

Sensitivity  (mM) 

0.99 

0.2 

0.62 

7 

Diameter  (nm) 

15 

10 

E.  Kjeang  et  al.  /  Journal  of  Power  Sources  158  (2006)  1-12 


5 


Fig.  2.  Normalized  Michaelis-Menten  curve  for  enzyme  kinetics. 


rates  are  then 


MDH  : 


Rt 


75 

~~  10"16  X  6.02  X  1023 
=  1.2  x  10-6  mol  m-2  s-1 


FalDH  : 


Ri 


35 

~~  10“16  x  6.02  x  1023 
=  5.7  x  10-7  molm-2  s_1 


rUH  :  /t  =  - tt - yv 

max  1(T16  x  6.02  x  1023 
=  1.9  x  10-7  mol  m-2  s_1 

If  all  the  electrons  released  by  the  chemical  reactions  are  cap¬ 
tured  by  the  electron  mediators,  the  maximum  current  density 
at  each  enzyme  patch  is  obtained  by 


operation  near  maximum  specific  activity  is  possible  even  with 
concentrations  almost  as  low  as  Km.  For  example  if  co  =  5Am, 
we  get  v  =  0.83t;max,  which  is  only  17%  short  of  the  maximum 
value. 

Specific  activities  found  in  the  literature  are  often  given  in 
terms  of  catalytic  turnover  numbers  per  unit  mass  of  protein 
in  the  enzyme  sample.  Assuming  that  the  enzymes  are  fully 
purified,  the  actual  molecular  turnover  number  can  be  calculated 
as 

^cat  =  oM  (5) 

The  enzyme  molecular  mass,  M,  is  also  available  in  enzyme 
databases.  Using  Eq.  (5)  the  maximum  turnover  numbers  may 
be  estimated  as  follows, 

A  10.4  x  1(T6  , 

MDH  :  4t,max  = - go -  X  430  x  106  =  75  s”1 

R  8.3  x  10“6  , 

FalDH  :  __  = -  x  250  x  106  =  35  s_1 

Cdi,mdx  zrrv 

r  2.43  x  10“6  , 

FDH  :  ^cat,max  = - go -  x  280  x  106  =  1 1  s”1 

We  denote  methanol,  formaldehyde  and  formate  oxidation  with 
superscripts  A,  B  and  C,  respectively.  Note  that  the  use  of  tab¬ 
ulated  values  for  these  calculations  can  be  treacherous,  as  the 
values  vary  depending  on  variations  in  molecular  structure  and 
ambient  conditions  such  as  temperature  and  pH.  Although  these 
turnover  numbers  are  expected  to  be  order-of-magnitude  cor¬ 
rect,  a  targeted  experimental  study  in  this  area  would  be  very 
beneficial.  To  determine  the  flux  of  substrate  and  product  to 
and  from  the  electrode  surface  for  a  reaction  rate  limited  sys¬ 
tem,  we  require  estimates  of  the  average  size  of  the  enzyme 
molecules  and  the  average  surface  coverage.  For  an  order-of- 
magnitude  estimate  we  assume  complete  monolayer  surface 
coverage  and  size  equivalent  to  the  cross-sectional  area  used 
for  ADH,  Ac  =  10-16m2,  for  all  three  enzymes.  The  turnover 


*max 


—  zFRmax 


(6) 


Here,  z  represents  the  number  of  electrons  released  by  oxidation 
of  one  substrate  molecule.  Note  that  all  turnover  parameters 
and  current  densities  also  depend  on  bulk  concentration  and 
enzyme  sensitivity  via  the  Michaelis-Menten  model  for  first- 
order  enzyme  catalysis  given  in  Eq.  (4).  The  actual  turnover 
rates  and  current  densities  are  thus  written  as 


R  =  R 


max 


CQ 

Am  +  Cp 


(7) 


i  =  i 


max 


CQ 

Am  +  CQ 


(8) 


As  before,  if  the  substrate  concentration  is  high  compared  to  the 
sensitivity  (co  »  Am),  the  Michaelis-Menten  kinetic  effects  can 
be  neglected.  This  is  not  always  the  case  in  consecutive  chemical 
processes,  as  some  concentrations  will  approach  zero,  especially 
if  effectively  complete  fuel  utilization  is  achieved. 

A  complete  biofuel  cell  must  also  contain  a  cathode,  where 
the  protons  and  electrons  released  at  the  anode  are  captured 
in  order  to  provide  a  closed  electrical  circuit.  The  need  for 
high  pressures  and  temperatures  make  conventional  02-reducing 
cathodes  incompatible  with  biocatalytic  anodes.  For  consistency 
with  the  specificity  of  the  biofuel  cell  concept,  the  cathodic 
chemical  reaction  should  also  be  catalyzed  by  suitable  redox 
enzymes.  Two  options  are  considered  in  this  work,  either 
reduction  of  oxygen  to  water  catalyzed  by  the  bilirubin  oxi¬ 
dase  enzyme  or  peroxide  reduction  to  water  catalyzed  by  the 
catalase-peroxidase  enzyme.  Peroxide  is  preferable  from  an  oxi¬ 
dant  supply  point-of-view  given  that  it  is  miscible  with  water 
in  all  proportions  and  can  easily  be  added  to  the  liquid  mix¬ 
ture  in  the  microchannel.  Oxygen  for  the  cathodic  reaction  can 
also  be  dissolved  in  water,  but  only  in  relatively  low  concen¬ 
trations.  Alternatively,  using  a  gas  permeable  backing  material 
such  as  poly(dimethylsiloxane)  (PDMS)  [27],  oxygen  from  the 
ambient  air  could  diffuse  through  the  backing  layer  to  supple¬ 
ment  the  dissolved  oxygen  for  the  cathodic  reaction.  Oxygen 
has  the  advantage  of  being  available  in  the  ambient  air  whereas 
peroxide  must  be  continuously  supplied  from  a  reservoir.  The 
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overall  cell  potential  is  however  limited  by  the  potential  dif¬ 
ference  between  the  anode  and  cathode  redox  couples,  and  it 
appears  that  the  theoretical  potential  difference  achievable  with 
peroxide  is  somewhat  higher  than  for  oxygen.  Moreover,  the  sin¬ 
gle  enzyme  cathodic  reactions  described  here  are  expected  to  be 
faster  than  the  anodic  reactions  [28,29].  The  focus  of  this  study 
is  the  biofuel  cell  anode,  and  its  inherent  kinetics  are  assumed 
to  be  independent  of  the  cathode  characteristics. 

2.2.  Species  transport 

The  diffusion  coefficient  of  small  molecules  in  water  at  room 
temperature  is  normally  around  D  ~  10-9  m2  s-1,  but  increases 
exponentially  with  temperature.  An  enzymatic  biofuel  cell  oper¬ 
ating  at  its  optimum  temperature,  say  40  °C,  is  likely  to  experi¬ 
ence  somewhat  faster  diffusion.  Kulikovsky  [30]  recommends 
a  reference  value  of  £>ref  =  1.58  x  10_9m2s_1  at  7j-ef  =  298K 
for  methanol  in  water  in  void  pores,  and  use  the  temperature 
dependence 

T>MeOH/H2o(T)  =  Aef  exp(0.026236(£  -  £ref))  (9) 

At  optimum  temperature  (40  °C)  we  get  £>  =  2.34  xlO-9 
m2  s-1.  This  value  is  assumed  to  be  valid  for  all  solutes. 

Total  species  flux  by  diffusion  and  advection  is  described  by 

J  =  —  £>Vc  +  cu  (10) 

The  ratio  of  advection  to  diffusion  is  defined  by  the  Peclet  num¬ 
ber 


The  average  bulk  velocity  is  denoted  by  u ,  and  £c  is  a  character¬ 
istic  length.  For  high  Peclet  numbers,  we  may  neglect  diffusive 
flux  at  the  inlet  and  outlet,  and  species  transport  becomes  propor¬ 
tional  to  the  local  concentration.  With  this  assumption,  a  useful 
analytical  expression  for  the  species  concentration  distribution 
may  be  derived.  Taking  heterogeneous  chemical  reactions  into 
account  and  assuming  maximum  reaction  rates,  the  species  bal¬ 
ance  in  a  channel  with  rectangular  cross-section  can  be  written 
as 

cmuhw  =  'y^R{nSiXUw  +  coutuhw  (12) 

j 

In  this  relationship,  cm  is  the  inlet  concentration,  u  the  average 
velocity,  h  the  channel  height,  w  the  channel  width,  £  the  length 
of  each  active  enzyme  patch  and  cou t  the  outlet  concentration. 
For  a  single  reaction  process  with  only  one  type  of  enzyme  fully 
covering  the  anode  surface,  the  species  balance  at  downstream 
position  x  in  the  microchannel  becomes 

C[nuh  =  £max-X  +  c(x)uh  (13) 

Solving  for  the  concentration,  we  get 

/  x  ^max  /1 

c(x)  =  cin - — —  x  (14) 

un 

This  analytical  result  will  be  used  to  assess  the  accuracy  of  the 
numerical  model. 


The  biofuel  cell  system  is  expected  to  exhibit  low  Damkohler 
numbers  (reaction  rate  limited)  and  therefore  essentially  con¬ 
stant  substrate  concentration  in  the  cross-stream  direction.  In 
order  to  get  complete  fuel  utilization,  we  can  then  design  the 
microchannel  to  match  zero  concentration  at  the  outlet.  But 
because  of  degrading  enzyme  kinetics  it  is  rather  inefficient  to 
use  concentrations  close  to  zero.  Moreover,  the  advective  flux 
is  reduced  for  near-zero  concentrations  and  the  diffusive  part 
should  again  be  added.  Thus  a  2D  computational  model  incor¬ 
porating  diffusive  flux  and  concentration  dependent  enzyme 
kinetics  is  required.  Fuel  utilization  based  on  total  flux  for  a 
single  reaction  process  is  calculated  as 

FU  =  1  -  —  (15) 

Tin 

To  determine  overall  fuel  utilization  for  a  consecutive  process 
we  need  to  include  the  other  reaction  steps  as  well.  For  methanol, 
two  electrons  are  released  in  each  step,  so  for  complete  fuel  uti¬ 
lization  we  get  six  electrons  per  methanol  molecule.  Overall  fuel 
utilization  in  terms  of  electrons  doing  actual  work  can  therefore 
be  defined  as 

FU  =  1 - ^  (16) 

6  J- 

u*yin 

For  instance,  100%  FU  for  the  first  step  in  a  three-step  process 
with  0%  FU  for  the  secondary  steps  would  result  in  an  overall 
FU  of  33%  using  this  equation. 

In  electrochemical  devices,  cross-stream  transport  may 
appear  due  to  water  involvement  in  the  chemical  reactions  and 
electro-osmotic  drag.  But  in  microscale  biofuel  cells  where  dif¬ 
fusion  is  fast  compared  to  reaction  rates  such  effects  do  not 
intervene  with  substrate  transport  and  may  readily  be  neglected. 

2.3.  CFD  model 

A  numerical  model  that  solves  the  Navier-Stokes  equations 
for  the  laminar  velocity  profile  of  the  bulk  flow  and  the  con¬ 
vective  diffusion  equation  for  solute  species  transport  is  imple¬ 
mented  in  the  commercially  available  software  package  FEM- 
LAB  3.1.  The  2D  channel  geometry  under  study  is  h  =  20  p,m  in 
height  and  L  =  30  mm  long,  which  is  considered  fixed  throughout 
the  analysis.  The  continuum  assumption  is  appropriate  for  liquid 
flows  in  microchannels  of  this  scale,  justifying  the  use  of  contin¬ 
uous  flow  theory  and  established  computational  fluid  dynamics 
methods  [31].  The  Navier-Stokes  equations  for  conservation  of 
mass  and  momentum  for  steady  incompressible  flow  are  given  as 

V  •  3  =  0  (17) 

pg-Vp  +  jiV2u  =  0  (18) 

No-slip  boundary  conditions  are  applied  at  the  channel  walls 
with  uniform  inlet  velocity  u  and  constant  outlet  pressure.  The 
flows  modeled  here  exhibit  low  Reynolds  numbers  (Re  ~  10-3), 
and  are  steady  and  fully  developed  over  the  majority  of  the 
computational  domain.  Gravity  effects  are  neglected,  and  the 
flow  is  assumed  to  be  isotropic  and  isothermal  throughout  the 
domain.  The  convective  diffusion  equation  for  solute  species 
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conservation,  with  passive  scalar  treatment  and  the  assumption 
of  infinite  dilution,  is  written  as 

dc* 

p - b  pu  •  Vc*  =  V(pDVc*)  (19) 

dt 

cl  represents  molar  concentration  of  solute,  and  one  separate 
convective  diffusion  equation  is  needed  for  each  of  the  partici¬ 
pating  species.  Source  terms  can  readily  be  ignored  since  there 
are  no  homogeneous  chemical  reactions.  Dirichlet  boundary 
conditions  are  used  at  the  inlet,  that  is,  constant  primary 
substrate  concentration  Cin  =  500mM  representing  a  fuel 
reservoir  and  zero  secondary  species  concentrations.  Neumann 
boundary  conditions  are  assigned  to  the  cathode  surface  at 
y  =  h,  describing  an  impermeable  wall  with  zero  concentration 
gradients  (dcl / dy)\y=h  =  0.  At  the  anode  surface  at  y  =  0  we 
have  concentration  dependent  flux  due  to  the  heterogeneous 
chemical  reactions.  Incorporating  Michaelis-Menten  kinetics 
for  the  enzyme  catalyzed  reactions,  the  boundary  conditions 
can  be  written  as  D(dcl /dy)\y=o  =  RlmSi^(cl / (Klm  +  c1))  for 
sinks  and  D(dcl  /  dy)\y=o  =  —  RLaxid /(K^  +  c-7))  for  sources 
with  i  /  j.  Purely  advective  flux  is  applied  for  all  solutes  at  the 
outlet,  (dcl /dx)\x=L  =  0. 

The  governing  equations  are  coupled  through  velocity  and 
chemical  reactions,  and  can  be  solved  either  sequentially  or 
simultaneously.  The  boundary  conditions  for  the  species  con¬ 
servation  equations  depend  on  concentration,  which  is  one  of 
the  outputs  of  each  iteration,  and  must  be  updated  accordingly. 
For  a  consecutive  process  the  anode  boundary  conditions  may 
also  depend  on  the  concentration  of  other  species.  For  example, 
the  second  species  in  a  two-step  process  is  produced  at  the  same 
rate  as  the  first  species  is  consumed,  i.e.  its  boundary  condition 
at  the  production  site  will  depend  on  the  reaction  kinetics  of 
the  first  species.  To  solve  all  the  governing  equations  of  a  con¬ 
secutive  system  simultaneously  thereby  becomes  very  memory 
intensive.  For  improved  efficiency,  it  is  preferable  to  solve  the 
equations  sequentially  as  follows: 

(1)  Navier-Stokes  equations  ->  u,  p 

(2)  u ,  cA  ->  convective  diffusion  eq.  for  species  A  — >  cA  (iter¬ 
ate) 

(3)  u ,  cA,  cB  — >  convective  diffusion  eq.  for  species  B  — >  cB 
(iterate) 

(4)  u ,  cB ,  cG  — >  convective  diffusion  equation  for  species 
C  — >  cG  (iterate) 

These  types  of  multiphysics  problems  are  conveniently  tack¬ 
led  using  a  finite  element  based  methodology,  which  can  handle  a 
wide  range  of  coupled  non-linear  physical  problems  with  depen¬ 
dent  boundary  conditions.  A  direct  solver  type  that  converts 
the  non-linear  governing  equations  to  linear  systems  is  applied, 
using  a  non-iterative,  fast  and  unconditionally  stable  (always 
converging)  approach. 

Mesh  optimization  is  an  important  part  of  problem  solv¬ 
ing  using  CFD.  Simple  elongated  unstructured  triangular  ele¬ 
ments  with  manually  refined  sections  are  used  for  this  study, 
mostly  because  of  their  flexibility.  A  thorough  mesh  convergence 


analysis  was  performed  for  simple  one-step  ethanol  oxidation 
catalyzed  by  ADH  enzymes.  The  substrate  concentration  was 
intentionally  kept  high  such  that  Michaelis-Menten  kinetics  can 
be  ignored.  Mesh  control  was  assessed  by  comparing  numer¬ 
ical  results  to  the  analytical  solution  derived  in  Section  2.2. 
The  selected  mesh,  including  local  mesh  refinements  in  regions 
with  high  gradients,  showed  an  accurate  concentration  profile 
as  compared  to  the  analytical  result  in  Eq.  (14)  and  was  also 
computationally  viable  in  the  context  of  direct  solver  memory 
limitations. 

3.  Results  and  discussion 

Maximum  current  densities  obtained  from  Eq.  (6)  are  /^ax  = 
4  pA  cm-2  for  FDH,  /B  ax  =  11  pA  cm-2  for  FalDH  and  /Aax  = 
24  pA  cm-2  for  MDH.  Ethanol  oxidation  by  ADH  can  reach 
a  maximum  current  density  of  =  50  pAcm-2,  which  is 
close  to  the  value  reported  by  Moore  et  al.  [14],  53  pA  cm-2  for 
Nafion  immobilized  enzymes.  However,  Moore’s  group  points 
out  that  their  biofuel  cell  anode  is  mass  transfer  limited  because 
of  slow  species  transfer  through  the  pore  network  of  the  Nafion 
structure,  and  does  not  fully  utilize  its  feasible  kinetics.  With 
electron  mediator  cross-linking  as  in  our  case,  enzymes  exposed 
directly  to  the  fuel  facilitate  much  faster  mass  transfer  and  con¬ 
sequently  higher  current  densities.  Using  enzyme  kinetics  data 
for  enzymes  in  solution,  as  in  this  model,  is  therefore  likely  to 
underestimate  reaction  rates,  and  provide  a  conservative  esti¬ 
mate  of  fuel  cell  performance. 

Fig.  3  presents  the  concentration  distribution  in  the  chan¬ 
nel  for  an  ethanol  fuelled  ADH  anode  with  the  selected  mesh, 
and  also  shows  the  extremely  large  aspect  ratio  of  the  geometry 
(1500:1).  The  inlet  velocity  is  in  this  case  tailored  for  50%  the¬ 
oretical  fuel  utilization.  The  substrate  concentration  decreases 
linearly  with  axial  length,  with  an  outlet  value  of  254  mM, 
near  the  predicted  250  mM  for  ~50%  fuel  utilization.  The  dis¬ 
crepancy  from  the  analytically  predicted  value  appears  when 
Michalis-Menten  kinetics  with  relatively  high  sensitivity  are 
introduced;  marginally  reducing  the  reaction  rates  even  at  high 
substrate  concentrations.  The  average  current  density,  deter¬ 
mined  by  boundary  integration,  is  approximately  equal  to  /AaxH. 
As  shown  by  the  contours  in  Fig.  3,  the  substrate  concentration 

c,  mM 


Fig.  3.  Ethanol  concentration  in  the  microchannel  with  an  ADH  catalyzed  anode 
at  y  =  0,  plotted  with  magnified  channel  height.  True  channel  geometry  is  shown 
above  (aspect  ratio  1500:1).  Conditions:  cm  =  500mM,  u  =  1.56  x  10-5  ms-1, 
h  =  20  x  10"6m,L  =  0.03  m. 
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is  roughly  constant  in  the  cross-stream  direction,  indicating  a 
reaction  rate  limited  process. 

To  estimate  the  required  pumping  power  for  the  microchannel 
flow,  we  calculate  the  pressure  drop  given  from  laminar  flow 
theory  [32]  as 


A  p  = 


32  iiLu 

D\ 


(20) 


With  hydraulic  diameter  =  2  h  the  pressure  drop  becomes 
Ap  =  5.9Pa.  The  pumping  power  is  equal  to  the  pressure  drop 
multiplied  by  the  volumetric  flow  rate 


Wp  =  V  A  p  =  uhwAp 


(21) 


Furthermore,  the  overall  power  output  of  the  fuel  cell  can  be 
estimated  as 


Wout  =  IVout  =  iLwVout  (22) 

Per  unit  width,  we  get  Wv/w  =  10-9Wm-1  and  Wout/w  = 
10-2  W  m_1  (Vout  =  1  V).  The  pumping  power  is  thereby  neg¬ 
ligible. 

A  similar  anode  model  is  also  established  for  three-step 
consecutive  methanol  oxidation  to  carbon  dioxide.  Six  useful 
electrons  are  released  in  this  process  as  opposed  to  only  two  for 
the  one-step  ethanol  process.  Considering  the  molecular  mass 
of  methanol  and  ethanol,  respectively,  the  energy  density  of 
methanol  is  much  higher.  In  fact,  4-5  times  more  useful  elec¬ 
trons  can  be  captured  per  unit  mass  of  fuel.  The  methanol  inlet 
concentration  is  set  to  500  mM  and  the  velocity  is  reduced  rel¬ 
ative  to  the  ADH  model  because  of  additional  steps  and  slower 
kinetics.  As  a  starting  point,  the  enzymes  are  distributed  along 
the  anode  surface  in  three  discrete,  equally  sized  zones  of  10  mm 
in  length  (case  1).  The  anode  surface  boundary  is  therefore 
divided  into  three  parts  and  the  turnover  rates  of  MDH,  FalDH 
and  FDH  are  implemented  in  a  stepwise  fashion  and  assumed 
to  be  completely  substrate  specific.  One  electrical  contact  is 
required  for  each  enzyme  zone  to  avoid  potential  losses,  and 
the  zones  have  to  be  sufficiently  separated  to  prevent  leak  cur¬ 
rents.  We  assume  this  spacing  to  be  thin  compared  to  the  zone 
lengths  and  apply  the  same  microchannel  geometry,  material 
properties  and  other  boundary  conditions  as  in  the  ADH  study. 
The  solid  lines  of  Fig.  4  present  the  resulting  concentration  dis¬ 
tributions  for  a  50%  methanol  fuel  utilization  predictive  inlet 
velocity.  This  figure  effectively  illustrates  the  reaction  rate  lim¬ 
ited  consecutive  anodic  process  in  a  microchannel.  The  methanol 
concentration  decreases  linearly  over  the  MDH  enzyme  zone, 
allowing  the  formaldehyde  product  concentration  to  increase 
at  the  same  rate.  Downstream  of  the  MDH  zone,  the  leftover 
methanol  concentration  remains  constant.  On  the  FalDH  enzyme 
surface,  formaldehyde  is  consumed  to  create  formate  on  a  one- 
to-one  basis.  The  slope  of  the  concentration  curves  is  less  steep, 
indicating  slower  kinetics  than  for  MDH.  Formate  is  slowly  con¬ 
sumed  on  the  FDH  surface  while  the  other  concentrations  are 
constant.  Carbon  dioxide  formed  on  this  surface  is  assumed  to  be 
dissolved  or  vented  off  in  a  way  that  does  not  disturb  the  kinetics 
of  the  FDH  enzymes.  The  substrate  concentrations  at  the  out¬ 
let  are  c^ut  =  277  mM,  c®ut  =  104  mM  and  c^ut  =  80  mM  for 


MDH  FalDH  FDH 


x,  mm 


Fig.  4.  Species  concentration  at  the  anode  surface  from  inlet  to  outlet.  Solid  lines 
indicate  results  for  an  inlet  at  x  =  0,  whereas  dashed  lines  correspond  to  results 
for  an  extended  inlet  section  (— 0.01  <x<0m).  Conditions:  cA  =  500 mM, 
c?  =  OmM,  cT  =  0  mM,  u  =  2.4  x  10-6  ms-1  ,h  =  20  x  10-6  m.  Enzyme  pat¬ 
terning:  La  =Lb  =LC  =  0.01  m  (case  1). 


the  initial  case  without  an  extended  inlet  section  (solid  lines). 
The  average  current  density  is  close  to  the  maximum  values  for 
the  individual  enzymes,  but  ideally,  the  outlet  concentrations 
should  all  be  zero.  Nonzero  quantities  indicate  unused  fuel  and 
incomplete  fuel  utilization,  which  in  this  case  is  FU  =  25%.  The 
reason  why  c^ut  is  much  larger  than  the  predicted  250  mM  is 
that  there  is  a  significant  diffusive  flux  (about  5%)  entering  the 
domain  at  the  inlet.  To  isolate  this  effect,  an  extended  1  cm  long 
non-catalyzed  section  is  added  to  the  inlet  region  in  the  numer¬ 
ical  model.  Results  with  an  extended  inlet  region  are  shown  as 
dashed  lines  in  Fig.  4  and  compared  to  the  previous  results.  The 
outlet  concentrations  are  now  c^ut  =  252  mM,  c^ut  =  128  mM 
and  c^ut  =  80  mM  with  overall  fuel  utilization  FU  =  27%.  The 
main  difference  is  that  all  species  are  allowed  to  diffuse  across 
the  interior  boundary  at  x  =  0,  even  though  the  methanol  flux 
entering  the  domain  at  v=—  0.001m  is  purely  advective.  In 
reality,  the  fuel  reservoir  and  electrodes  ought  to  be  separated 


MDH  FalDH  FDH 


x,  mm 


Fig.  5.  Species  concentration  at  the  anode  surface  from  inlet  to  outlet  for 
various  average  velocities  (u  =  2.4  x  10-6ms-1  for  dotted  lines,  u  =  0.8  x 
10-6ms-1  for  dashed  lines  and  u  =  0.3  x  10-6ms-1  for  solid  lines).  Enzyme 
patterning:  LA  =  3.2  mm,  LB  =  6.7  mm,  Lc  =  20. 1  mm  (case  2). 
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Fig.  7.  Concentration  distributions  at  the  anode  surface  for  two  bulk  velocities; 
high(n  =  8x  10-7  ms-1)  and  low  (u  =  3  x  10-7  ms-1).  Enzyme  patterning: 
mixed  1:1:1  (case  3). 

in  this  manner  by  some  inlet  feed  section  with  no  chemical 
reactions. 

In  order  to  improve  fuel  utilization,  we  have  to  take  advan¬ 
tage  of  the  difference  in  turnover  rates  for  the  various  enzymes. 
The  zones  with  low  turnover  rates  should  be  extended  at  the 
expense  of  the  high  rate  zones.  Since  each  reaction  generates  an 


equivalent  number  of  electrons,  the  optimal  lengths  of  the  zones 
correspond  to  the  ratio  l/RA :  l/RB :  l/Rc .  Assuming  that  R  =  RmSiX 
everywhere,  we  get 

V  =  (7?max)  L,  i  =  A,B,C,  j  =  A,  B,  C  (23) 

y^j(^max) 

For  the  three- step  methanol  process,  the  optimal  lengths  become 
La  =  3.2  mm,  LB  =  6.7  mm  and  Lc  =  20. 1  mm,  which  are  applied 
in  case  2.  Fig.  5  shows  concentration  profiles  with  respect  to 
these  zone  lengths,  obtained  for  three  different  bulk  velocities. 
Again,  an  extended  inlet  section  is  used.  For  the  same  velocity  as 
in  case  1  (dotted  lines),  the  methanol  concentration  remains  high 
and  the  other  concentrations  are  low  throughout.  Case  1  does 
actually  perform  better  than  case  2  at  high  velocities.  In  order 
to  reach  high  fuel  utilization  the  velocity  must  be  reduced,  as 
can  be  seen  by  the  improved  methanol  concentration  drop  as  the 
velocity  is  cut  down.  Simultaneously,  as  we  reduce  the  velocity, 
Pe  numbers  approach  unity  and  an  increased  fraction  of  diffu¬ 
sive  flux  enters  the  non-catalyzed  inlet  section.  An  important 
issue  to  note  in  the  case  of  low  Pe  is  that  the  secondary  species 
are  able  to  diffuse  back  into  the  inlet  section.  For  the  lowest 
velocity  (u  =  0.3  x  10-6ms-1),  both  formaldehyde  and  for¬ 
mate  are  able  to  diffuse  back  all  the  way  to  the  fuel  reservoir. 
Back  diffusion  of  secondary  species  will  in  principle  have  no 
effect  on  fuel  utilization,  but  it  can  reduce  current  densities  due 


Fig.  8.  Current  densities  (a)  and  fuel  utilization  (b)  vs.  Peclet  number.  Enzyme  patterning:  mixed  1:1:1  (case  3). 
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to  a  lack  of  secondary  fuel  at  the  enzyme  sites.  Current  densities 
at  the  individual  enzyme  zones  as  well  as  the  overall  current 
density  are  shown  in  Fig.  6a,  where  the  Peclet  numbers  are 
based  on  the  characteristic  length  scale  for  streamwise  diffu¬ 
sion,  Lc  =  0.01  m.  The  current  densities  decrease  rapidly  as  the 
rate  of  diffusion  overrides  the  rate  of  advection  (Pe  <  1).  On 
the  other  hand,  all  current  densities  are  remarkably  stable  in  the 
advection  dominated  regime.  Fuel  utilization  is  near  100%  in  the 
diffusion  dominated  and  mixed  regimes,  as  indicated  in  Fig.  6b, 
unaffected  by  back  diffusion.  It  is  also  shown  that  the  enzyme 
turnover  rates  cannot  keep  up  with  the  species  transport  as  the 
velocity  reaches  some  threshold  value,  causing  less  fuel  efficient 
operation  with  unused  fuel  at  the  outlet. 

Another  option  for  the  anode  enzyme  patterning  is  to  mix 
all  three  types  throughout  the  entire  electrode  surface  using 
only  one  electrical  contact,  which  is  preferable  from  a  fabri¬ 
cation  point-of-view.  The  simplest  option  is  to  mix  the  enzymes 
on  a  1:1:1  molecular  basis,  denoted  by  case  3,  such  that  each 
enzyme  covers  one  third  of  the  surface  area  and  the  turnover  rates 
become 


max 


C 

A  m  +  C 


(24) 


Concentration  profiles  for  case  3  are  presented  in  Fig.  7  for 
two  different  bulk  velocities.  In  general,  the  concentration  distri¬ 
butions  are  similar  to  the  results  of  case  1 ,  where  a  large  flux  of 
secondary  species  is  leaving  the  domain  at  the  outlet.  To  measure 
this  effect,  current  densities  (a)  and  overall  fuel  utilization  (b)  are 
given  in  Fig.  8  as  functions  of  the  Peclet  number.  High  current 
densities  are  obtained  in  the  advection  dominated  regime  but  the 
fuel  utilization  for  high  Pe  is  poor.  The  fuel  efficiency  is  only  suf¬ 
ficient  in  the  diffusion  dominated  regime.  It  is  therefore  difficult 
to  find  an  appropriate  balance  between  these  two  parameters 
for  the  mixed  enzyme  electrode  design,  unless  the  individual 
enzyme  surface  coverage  ratios  are  modified.  Matching  the  sur¬ 
face  coverage  to  the  turnover  rates,  as  in  case  2,  we  utilize  the 
ratios  l/RA:l/RB:l/Rc.  The  new  turnover  rates  adopted  in  case 
4  are  calculated  by 


R  = 


(^max) 


-1 


£,-(*  i 


max) 


_  j  ^max 

Krr 


+  c 


j  =  A,  B,  C 


(25) 


Fig.  9.  Concentration  distributions  at  the  anode  surface  for  two  bulk  velocities; 
high  (u  =  8  x  10-7  ms-1)  and  low  (u  =  3  x  10-7  ms-1).  Enzyme  patterning: 
mixed  l/RA:l/RB:l/Rc  (case  4). 


The  concentration  profiles  obtained  for  case  4  are  shown  in  Fig.  9 
for  the  same  velocities  as  in  Fig.  7.  The  rate  of  methanol  con¬ 
sumption  is  slower  than  in  case  3  as  the  MDH  fraction  is  reduced, 
but  significant  improvements  are  attained  in  terms  of  secondary 
species  outlet  concentrations.  Fig.  10  presents  current  densi¬ 
ties  and  fuel  utilization  for  case  4.  In  this  case,  fuel  utilization 
remains  high  and  steady  for  a  wider  range  of  Peclet  numbers. 
Maximum  current  densities  are  lower  than  for  case  3  but  can  be 
realized  for  most  Pe  above  unity.  It  is  possible  to  combine  high 
current  densities  with  high  fuel  usage  for  this  enzyme  patterning 
strategy.  For  example,  Pe  =  1.5  provides  7.0  pAcm-2,  which  is 
99%  of  the  maximum  current  density  for  case  4,  at  99%  fuel 
utilization. 

To  summarize  the  results  obtained  for  the  anodic  three- step 
methanol  oxidation  process,  Fig.  1 1  presents  a  comparison  of 
current  densities  and  fuel  utilization  for  cases  1-4  at  low  Peclet 
numbers.  The  diffusion  dominated  regime  ( Pe<  1)  is  not  rec¬ 
ommended  due  to  low  current  densities.  Fuel  cell  performance 
is  generally  better  in  the  mixed  regime  (1  <Pe<2)  where  rea¬ 
sonably  high  current  densities  are  met  by  nearly  complete  fuel 
utilization,  with  case  4  showing  superior  characteristics.  In  the 


Fig.  10.  Current  densities  (a)  and  fuel  utilization  (b)  vs.  Peclet  number.  Enzyme  patterning:  mixed  l/RA:l/RB:l/Rc  (case  4). 
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Fig.  11.  Current  density  (a)  and  fuel  utilization  (b)  comparison  for  enzyme  arrangements  according  to  cases  1-4. 


advection  dominated  regime,  cases  1  and  3  are  the  preferable 
choices  although  fuel  utilization  is  poor. 

4.  Conclusions 

A  micro  structured  enzymatic  biofuel  cell  architecture  for 
consecutive  chemical  reactions  was  proposed,  and  species  trans¬ 
port  coupled  with  laminar  flow  and  Michaelis-Menten  kinetics 
was  studied  in  2D  using  a  numerical  approach.  Biofuel  cell 
performance  was  shown  to  be  limited  by  the  reaction  rates 
associated  with  enzyme  kinetics.  Turnover  rates  for  individ¬ 
ual  enzymes  were  key  parameters  throughout  this  analysis  and 
directly  determined  the  realizable  current  densities.  The  pump¬ 
ing  power  required  for  the  microchannel  flow  was  determined 
to  be  negligible  compared  to  the  predicted  output  power  of  the 
unit  cell. 

The  current  densities  obtained  for  one-step  enzyme  catalyzed 
ethanol  oxidation  were  generally  higher  than  for  consecutive 
methanol  oxidation.  If  maximum  current  density  is  the  only 
target,  one-step  ethanol  oxidation  is  thereby  recommended. 
Methanol  is,  however,  the  more  attractive  fuel  in  terms  of  energy 
density,  as  4-5  times  more  electrons  can  be  captured  per  unit 
mass  of  fuel.  Four  separated  and  mixed  electrode  enzyme  pat¬ 
terning  strategies  were  presented,  and  tested  with  various  bulk 
velocities  to  optimize  overall  current  density  and  fuel  utiliza¬ 
tion.  Nearly  complete  fuel  utilization  was  accomplished  in  the 
diffusion  dominated  and  mixed  species  transport  regimes.  Local 
current  densities  are  high  and  steady  in  the  advection  dominated 
regime,  but  become  significantly  reduced  for  low  velocities 
because  of  back  diffusion  and  fuel  supply  shortcomings.  The 
mixed  transport  regime  is  therefore  particularly  attractive  for 
biofuel  cell  operation,  with  particularly  superior  characteristics 
shown  for  mixed  enzyme  patterning  that  accounts  for  individual 
turnover  rates.  Overall  current  density  is  mainly  limited  by  the 
slow  turnover  rates  for  formate  and  formaldehyde  oxidation  by 
FalDH  and  FDH.  Finding  new  or  modified  enzymes  with  higher 
turnover  numbers  is  a  sound  strategy  to  improve  the  perfor¬ 
mance  of  this  technology.  Inclusion  of  proper  enzyme-electrode 
tethering  is  expected  to  enhance  current  densities,  and  accu¬ 
rate  measurements  of  individual  turnover  numbers  for  tethered 
enzymes  are  desired  for  further  advancements. 
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